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INTRODUCTION 


la  recent  years  the  Interest  In  plastic  deformation  has  Increased 
considerably  In  research  as  well  as  in  practical  applications.  The  Importance 
of  realistic  constitutive  equations  has  been  emphasized  in  several  research 
workshops  (refs  1-3),  and  many  computer  codes  have  been  developed  as  seen  in  a 
recent  survey  paper  (ref  4).  According  to  time-independent  plasticity  theory, 
the  response  of  a  strain-hardening  material  is  specified  by  an  initial  yield 
condition,  a  hardening  rule,  and  a  flow  rule.  While  there  is  general  agree¬ 
ment  in  the  literature  over  which  initial  yield  condition  and  flow  rule  should 
be  used,  there  is  no  such  accord  with  regard  to  the  hardening  rule.  Most 
hardening  rules  in  present  use  are  well  documented  and  reviewed  in  the 
literature  (ref  5).  The  stress-strain  curves  in  Figures  l  and  2  illustrate 
the  fact  that  all  of  the  theories  are  capable  of  treating  the  monotonic 
loading  situation.  For  reversed  loading,  it  has  been  concluded  that  kinematic 
and  isotropic  hardening  models  represent  the  limits  of  the  actual  behavior, 
whereas  the  remaining  theories  are  capable  of  falling  anywhere  within  these 
limits  (refs  5,6).  This  is  certainly  true  for  many  strain-hardening  materials 
as  shown  in  Figure  1,  but  not  true  for  some  small  strain-hardening  materials 
as  shown  in  Figure  2.  The  experimental  works  by  Milligan,  Koo,  and  Davidson 
(ref  7)  and  Swift  (ref  8)  demonstrate  that  the  stress-strain  curves  for  a  high 
strength  steel  and  six  spring  steels  behave  like  that  shown  in  Figure  2.  None 
of  the  existing  theories  can  represent  this  material  behavior  reasonably  well 
for  two  reasons.  First,  the  initiation  of  reverse  yielding  occurs  much 

References  are  listed  at  the  end  of  this  report. 
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earlier  then  predicted  by  the  theories.  Second,  the  assumption  of  same  slope 
for  forward  and  reversed  loading  Is  not  valid. 

In  this  report  a  theoretical  model  Is  proposed  with  an  attempt  to  give  a 
better  representation  of  the  actual  material  behavior  especially  that  shown  In 
Figure  2  (refs  7,8).  The  Bauschinger  effect  factor  Is  treated  as  a  function 
of  overstrain.  The  strain-hardening  effect  la  taken  into  account  with 
different  parameters  used  for  forward  and  reversed  loading  processes.  The 
application  of  this  model  to  the  torsion  problem  in  a  cylindrical  bar  is 
reported. 

THEORETICAL  MODEL 

The  stress-strain  curve  during  loading  for  a  small  strain-hardening 
material  can  be  replaced  with  sufficient  accuracy  by  a  bilinear  elastic- 
plastic  model  as  shown  in  Figure  3.  For  the  plastic  portion,  the  yield  shear 
stress  t  is  related  to  the  plastic  shear  strain  yP  by 

t/t o  *  l  +  m?/(l-m)  and  5  -  (G/t0)yP  (1) 

where  is  the  initial  yield  shear  stress,  G  is  the  shear  modulus,  and  mG  is 
the  slope  of  shear  stress-strain  curve  in  the  plastic  portion. 

In  most  of  the  plasticity  theories,  the  curve  of  reverse  loading  is 
uniquely  defined  by  the  curve  of  the  first  loading.  The  present  model  does 
not  assume  such  a  relationship.  The  experimental  stress-strain  curve  during 
reversed  loading  will  be  used  directly.  A  piecewise  linear  representation  can 
be  used,  but  only  a  bilinear  approximation  is  chosen  here  as  shown  in  Figure 
3.  Choosing  a  new  coordinate  system  (t* ,  y*  )  with  origin  at  the  point  before 


unloading,  we  have  for  the  plastic  portion  of  the  reverse  yielding  curve 
t’/t0  -  t0'/t0  +  n'c7<l-m')  and  C  -  (G/t0)Y'P 
where  t0'  is  the  linear  drop  In  shear  stress  until  reverse  yielding  begins, 
m'G  is  the  slope  of  the  reverse  yielding  curve,  and  y'P  is  the  additional 
plastic  shear  strain  during  reversed  loading.  Taking  into  account  strain- 
hardening,  the  definition  of  the  Bauschinger  effect  factor  (BE?)  is 


BEF  -  (V-Ti)/*!  -  fUx)  (3) 

where  Tj  is  the  shear  stress  Just  before  unloading  occurs.  The  Bauschinger 
effect  factor  is  a  function  of  prestrain  as  shown  in  Figure  4.  This 
information  is  available  from  the  experimental  data  for  a  high  strength  steel 
(ref  7). 

According  to  equations  (1)  and  (3),  t0’  can  be  expressed  as  a  function  of 
plastic  strain  (C^)  just  prior  to  unloading  by 

I0'/T0  «  U  +  «C/(l-B)J[l  +  f(C)J  -  g(C)  (4) 

The  present  model  Is  quite  general  because  different  parameters  (m,m’ )  and 
functions  f(C)  can  be  used.  It  seems  interesting  to  discuss  two  special  cases 
of  the  present  material  model: 

(a)  m  ■  m' ,  f  -  1:  reduces  to  isotropic  hardening  model. 

(b)  m  ■  m* ,  g  ■  2:  reduces  to  klnematical  hardening  model. 

TORSION  BARS 

For  a  circular  bar  under  torsional  loading,  the  external  torque  M  can  be 
computed  by 


1 

M  ■  M/Mo  -  4  /  (T/T0)52d6 
0 


(5) 
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where  Hq  ■  t/2  «9t0,  (  -  r/a,  t  la  the  current  shear- • trees  corresponding  to 
location  r,  end  e  Is  the  radius  of  the  cross-section.  Geometric 
considerations  show  that  radial  lines  have  to  remain  straight  after 
deformation.  Thus,  one  concludes  that 

Y  -  YaC  -  or  (6) 

where  Ya  is  the  strain  at  the  outermost  fiber  and  a  is  the  angle  of  twist  per 
unit  length.  Since  a  yield  stress  is  introduced,  an  elastic  core  always 
exists  during  deformation  whose  radius  c  is  given  by 

c/a  -  Y0/Ya  and  Yc  -  t0/G  (7) 

If  a  material  is  linear  strain  hardening,  a  closed-form  loading  solution 
the  elastic-plastic  ractge  can  be  obtained  (ref  9).  The  plastic  shear  strain 
and  shear  stress  are  given  by 

5  -  <l-m)(Y/Y0-l>  (8) 

t/t0  -  (l-m)  +  mY/Y0  (9) 

Using  equations  (5)  through  (9),  we  can  compute  the  torque.  The  result  is 

-  me  l-m  , ,  , 

M  ■  —  + - (4-c3/a3)  (10) 

c  3 

Note  that  M  «  1  for  c  *  a  (initial  yielding),  and  M  »  (l/3)(4-c3/a3)  for 
perfectly  plastic  material. 

REVERSED  LOADING 

If  we  remove  the  external  torque  after  the  bar  is  twisted  beyond  the 
elastic  limit,  there  is  a  residual  angle  of  twist  (a").  The  ratio  of  the 
residual  angle  to  the  original  angle  of  twist  (a)  is  defined  as  the  spring- 
back  ratio  (a"/a).  Assuming  complete  elastic  unloading  (ref  9),  we  can 


r 


determine  o’,  end  the  spring-back  ratio  (o-oi')/a  is  given  by 


4  Y0  l  Y0 

a“/a  -  (l-m)[l - +-(--)"]  (11) 

3  Ya  3  ya 


This  formula  for  spring-back  ratio  can  be  used  if  reverse  yielding  will  not 
occur.  It  can  be  shown  that  this  is  true  when  the  following  inequality  Is 
satisfied: 


ma  1-m 


c 
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—  +  (4  --•;)<  g(  Ca)  (12) 

c  3  a* 

where  Ca  is  the  dimensionless  plastic  shear  strain  at  the  outermost  fiber  just 
prior  to  unloading.  The  upper  limit  for  a  or  Y  and  the  corresponding  M  can  be 
found  from  the  above  inequality. 

As  an  example,  when  m  *  0.1,  c/a  ■  0.1,  g  ■  2  (kinematic  hardening), 
inequality  (12)  is  not  satisfied,  l.e.,  reverse  yielding  occurs  during  the 
removal  of  the  external  torque  M  ■  2.1996 

The  occurrence  of  reverse  yielding  depends  on  the  initial  loading, 
reversed  loading,  and  the  material  models  Even  if  reverse  yielding  may  not 


occur  during  the  removal  of  external  torque,  it  may  still  occur  during  the 
reversed  loading.  Let  a'  be  the  unit  angle  of  twist  in  the  reversed  direction 
during  reversed  loading.  Geometric  considerations  again  lead  to 

V  -  Y.'C  ■  o'r  (13) 

where  Ya*  is  the  shear  strain  at  the  outermost  fiber  during  reversed 
loading.  Let  d  be  the  elastic-plastic  boundary.  In  the  elastic  tone 
(r  <  d),  we  have  Y*P  -  0  and  t*  -  GY'.  In  the  plastic  tone  (r  <  d),  the 
shear  stress  and  plastic  shear  strain  during  reverse  loading  can  be 
computed  by 


t'/t0  -  (l-m')g(C)  +  m’Y'/Yo 


(14) 
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and 


O  -  (l-m’JlY'/Yo  -  g ( O]  (16) 

AC  the  elastic-plastic  boundary  d,  *  0  and  Cd  =*  (1-ra) ( Yd/Y0~l)  which  lead 
to 

Yd'/Yo  =  g(«d>  (16) 

«d  =  (l-m)[(d/a)(Ya/Y0)  -  l]  (17) 

Using  equations  (13),  (16),  and  (17),  we  obtain 

d  gUd>  1  +  Cd/(l“m) 

-  - - = -  (18) 

a  Ya'/Y0  Ya/Y0 


Therefore,  we  can  compute  Yj  and  d  for  given  values  of  Ya  and  Ya' •  The  ratio 
Ya'/Ya  Is  related  to  4 d> 

Ya' / Ya  -  a'/a  -  g(5d)/[l  +  Cd/U"®)]  (19) 

Now  we  can  calculate  the  external  torque  M*  during  reversed  loading  by 

M*  -  M«/M0  -  (Ya’/Yo)!*'  +  (!-■’ )(d/a) **)  +  (1-m1)  j]  g(C)C2d S  (20) 


If 


KO  “  f0  +  El^  +  •••  +  fn-l  f*n_1 


then 

g(0  -  g0  +  glC  +  •••  +  SnO 


(21) 


Since  C  Is  linear  In  C,  we  can  carry  out  the  integration  in  equation  (20) 
explicitly  In  terms  of  d/a,  ra,  and  Ya/Y0«  In  general,  H'  is  a  function  of  Ya 
and  Ya'  with  parameters  a,  0,  xQ,  ra,  ra* ,  f0,  The  dependence  of  M' 

on  d/a  has  been  eliminated  through  equations  (18)  and  (19). 
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NUMERICAL  RESULTS  AND  DISCUSSIONS 

Consider  a  circular  bar  which  was  twisted  to  reach  Ya  »  10yo  where  Y0  * 
t0/G.  The  elastic-plastic  boundary  was  calculated  by  equation  (7)  to  be  c/a  * 
0.1.  The  applied  external  torque  can  be  calculated  by  equation  (10);  and  the 
results  are  M/Mq  -  1.333,  1.41967,  2.1996  for  ra  =*  0.0,  0.01,  0.1, 
respectively,  and  Mq  «  x/2  a3t0. 

During  the  removal  of  the  external  torque,  the  occurrence  of  reverse 
yielding  still  depends  on  the  material  model.  If  the  kinematic  hardening 
model  Is  used  with  m  =  m'  *  0.1,  then  reverse  yielding  will  occur  because 
Inequality  (12)  is  violated.  By  equating  M*  to  M  as  given  by  equations  (10) 
and  (20),  and  solving  for  the  elastic-plastic  boundary,  we  have  d/a  »  0.3910. 
The  residual  shear  strain  is  Ya'  *  2.2446yq  and  the  spring-back  ratio  during 
elastic-plastic  unloading  is  0.77554,  According  to  equation  (11)  on  the 
assumption  of  elastic  unloading,  the  spring-back  ratio  would  be  0.78003.  If 
ra  =  in*  *  0.01  is  used  with  either  kinematic  (g  =  2)  or  lostropic  (f  =»  1) 
hardening  models,  reverse  yielding  will  not  occur  during  unloading.  However, 
reverse  yielding  may  still  occur  If  the  Bauschlnger  effect  factor  is  small, 
e.g.,  f  »  0.3. 

When  the  external  torque  is  removed  and  then  applied  In  the  reversed 
direction,  reverse  yielding  will  occur  no  matter  what  material  models  are 
used.  Different  material  models  predict  quite  different  results.  Let  us 
consider  only  the  ^ase  of  Ya  -  Ya'  *  10yo  with  m  -  0.01.  The  external 
torque  to  reach  the  Initiation  of  reverse  yielding  will  be  -0.76033,  -0.58033, 
-0.24633  Mq  according  to  Isotropic  hardening,  klnemsttc  hardening,  the  present 
model  with  f  “  0.4,  respectively.  Further  reversed  loading  depends  also  on 
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the  slope  of  strain-hardening.  According  to  either  the  Isotropic  or  the 
kinematic  hardening  model,  the  slope  during  reverse  loading  is  the  same  as 
during  loading,  i.e.,  m'  ■  m  ■  0.01.  The  present  model  does  not  make  such  an 
assumption.  Different  slopes  such  as  m'  »  0.1,  0.3  can  be  used,  and  f(C)  can 
be  a  general  function  of  prestrain.  To  reach  Ya'  ■  Ya  ■  10Yo  (twist  the  bar 
forward  and  backward  in  the  same  amount),  we  need  to  apply  the  external  torque 
In  the  reversed  direction.  The  predicted  results  are  -1.31504,  -1.48844  Mq 
according  to  the  kinematic  hardening  (g  ■  2)  and  Isotropic  hardening  (f  -  1) 
models,  respectively.  If  the  Bauschinger  effect  factor  f  ■  0.4  is  used  with 
m*  *  0.01,  0.1,  0.3,  then  the  external  torque  will  be  -0.64716,  -1.3686, 
-2.97104  Mq,  respectively.  Therefore,  a  high  strength  steel,  with  m  ■  0.01,  f 
*  0.4,  and  ra'  -  0.3  requires  a  much  larger  torque  in  the  reversed  direction 
(-2.97104  Mo)  than  that  predicted  by  the  other  models,  even  though  the 
initiation  of  reverse  yielding  occurs  at  a  much  smaller  torque  (-0.24633  Mq) 
than  the  other  predictions. 

Since  different  material  models  predict  quite  different  results,  reverse 
loading  experiments  are  needed  for  the  characterization  of  the  material 
behavior  and  for  the  comparison  with  the  theoretical  predictions. 

SUMMARY  AND  RECOMMENDATIONS 

1.  The  proposed  model  for  reverse  yielding  can  give  a  better 
representation  of  the  actual  stress-strain  curve.  All  popular  plasticity 
models  fail  for  reverse  yielding  in  a  high  strength  steel  or  spring  steel. 

2.  The  present  model  has  been  applied  to  the  torsion  problem  in  a 
cylindrical  bar.  Por  a  high  strength  steel  bar,  the  results  indicate  that  the 


1 


initiation  of  reverse  yielding  occurs  at  a  much  smaller  torque  in  the  reversed 
direction  than  that  predicted  by  the  other  models.  Further  reversed  loading 
requires  a  much  larger  rate  of  Increase  In  the  applied  torque  than  the  other 
mode is . 

3.  More  research  works  to  develop  a  better  plasticity  theory  are  needed 
for  reversed  loading  and  cyclic  loading  problems.  Basic  experiments  are 
needed  to  characterize  the  material  behavior  and  to  compare  with  the 
theoretical  predictions. 
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CHIEF,  ENGINEERING  SUPPORT  BRANCH 

ATTN:  DRSMC-LCB-S  1 

-SE  I 

CHIEF,  RESEARCH  BRANCH 

ATTN:  DRSMC-LCB-R  2 

-R  (ELLEN  FOGARTY)  1 

-RA  1 

-RK  2 

-RP  1 

-RT  l 

TECHNICAL  LIBRARY  5 

ATTN:  DRSMG-LCB-TL 

TECHNICAL  PUBLICATIONS  &  EDITING  UNIT  2 

ATTN:  QRSMC-LCB-TL 

DIRECTOR,  OPERATIONS  DIRECTORATE  1 

DIRECTOR,  PROCUREMENT  DIRECTORATE  1 

DIRECTOR,  PRODUCT  ASSURANCE  DIRECTORATE  l 


NOTE:  PLEASE  NOTIFY  DIRECTOR,  BENET  WEAPONS  LABORATORY,  ATTN:  DRSMC-LCB-TL , 
07  ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST 


NO.  OF  NO.  OF 

COPIES  COPIES 


ASST  SEC  OF  THE  ARMY 

RESEARCH  &  DEVELOPMENT 

ATTN:  DEP  FOR  SCI  &  TECH  1 

THE  PENTAGON 

WASHINGTON,  D.C.  20315 

COMMANDER 

DEFENSE  TECHNICAL  INFO  CENTER 
ATTN:  DTIC-DDA  12 

CAMERON  STATION 
ALEXANDRIA,  VA  22314 

COMMANDER 

US  ARMY  MAT  DEV  &  READ  COMD 
ATTN:  DRCDE-SG  1 

5001  EISENHOWER  AVE 
ALEXANDRIA,  VA  22333 


COMMANDER 

ARMAMENT  RES  &  DEV  CTR 
US  ARMY  AMCCOM 

ATTN:  DRSMC-LC(D)  1 

DRSMC-LCE(D)  1 

DRSMC-LCM(D)  (BLDG  321)  1 

DRSMC-LCS(D)  1 

DRSMC-LCU(D)  1 

DRSMC-LCW(D)  1 

DRSMC-SCM-0  (PLASTICS  TECH  1 

EVAL  CTR, 

BLDG.  351N) 

DRSMC-TSS(D)  (STINFO)  2 

DOVER,  NJ  07801 


DIRECTOR 

BALLISTICS  RESEARCH  LABORATORY 
ARMAMENT  RESEARCH  &  DEV  CTR 
US  ARMY  AMCCOM  1 

ATTN:  DRSMC-TSB-S  (STINFO) 

ABERDEEN  PROVING  GROUND,  MD  21005 

MATERIEL  SYSTEMS  ANALYSIS  ACTV 

ATTN:  DRSXY-MP  1 

ABERDEEN  PROVING  GROUND,  MD  21005 


COMMANDER 
US  ARMY  AMCCOM 

ATTN:  DRSMC-LEP-L( R)  1 

ROCK  ISLAND,  IL  61299 

COMMANDER 

ROCK  ISLAND  ARSENAL 

ATTN:  SMCRI-ENM  (MAT  SCI  DIV)  1 

ROCK  ISLAND,  IL  61299 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENG  ACTV 
ATTN:  DRXIB-M  1 

ROCK  ISLAND,  IL  61299 

COMMANDER 

US  ARMY  TANK-AUTMV  R&D  COMD  1 

ATTN:  TECH  LIB  -  DRSTA-TSL 

WARREN,  MI  48090 

COMMANDER 

US  ARMY  TANK-AUTMV  COMD  1 

ATTN:  DRSTA-RC 

WARREN,  MI  48090 

COMMANDER 

US  MILITARY  ACADEMY 

ATTN:  CHMN,  MECH  ENGR  DEPT  1 

WEST  POINT,  NY  10996 

US  ARMY  MISSILE  COMD 
REDSTONE  SCIENTIFIC  INFO  CTR  2 

ATTN:  DOCUMENTS  SECT,  BLDG.  4484 

REDSTONE  ARSENAL,  AL  35898 

COMMANDER 

US  ARMY  FGN  SCIENCE  &  TECH  CTR 
ATTN:  DRXST-SD  1 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER, 
US  ARMY  AMCCOM,  ATTN:  BENET  WEAPONS  LABORATORY,  DRSMC-LCB-TL , 
WATERVLIET,  NY  12189,  OF  ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST  (CONT'D) 


NO.  OF 

NO.  OF 

COPIES 

COPIES 

COMMANDER 

DIRECTOR 

US  ARMY  MATERIALS  &  MECHANICS 

US  NAVAL  RESEARCH  LAB 

RESEARCH  CENTER 

2 

ATTN:  DIR,  MECH  DIV 

1 

ATTN:  TECH  LIB  -  DRXMR-PL 

CODE  26-27,  (DOC  LIB) 

1 

WATERTOWN,  MA  01272 

WASHINGTON,  D.C.  20375 

COMMANDER 

COMMANDER 

US  ARMY  RESEARCH  OFFICE 

AIR  FORCE  ARMAMENT  LABORATORY 

ATTN:  CHIEF,  IPO 

1 

ATTN:  AFATL/DLJ 

1 

P.O.  BOX  12211 

AFATL/DLJG 

1 

RESEARCH  TRIANGLE  PARK,  NC  27709 

EGLIN  AFB,  FL  32542 

COMMANDER 

METALS  &  CERAMICS  INFO  CTR 

US  ARMY  HARRY  DIAMOND  LAB 

BATTELLE  COLUMBUS  LAB 

1 

ATTN:  TECH  LIB 

t 

505  KING  AVENUE 

2800  POWDER  MILL  ROAD 

COLUMBUS,  OH  43201 

ADELPHIA,  MD  20783 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CTR 
ATTN:  TECHNICAL  LIBRARY  1 

CODE  X212 

OAHLGREN,  VA  2 2448 


NOTE:  PLEASE  NOTlrY  COMMANDER,  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER, 
US  ARMY  AMCCOM,  ATTN:  BENET  WEAPONS  LABORATORY,  DRSMC-LCB-TL, 
WATERVLIET,  NY  12189,  OF  ANY  ADDRESS  CHANGES. 


